Electrical detection of coherent 31 P spin quantum states 
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In recent years, a variety of solid-state qubits has been realized, including quantum dots 
superconducting tunnel junctions 0-01 and point defects 0.0. Due to its potential compatibility 
with existing microelectronics, t he p roposal by Kane 0,0 based on phosphorus donors in Si has also 
been pursued intensively |SL llOL lllj - A key issue of this concept is the readout of the 31 P quantum 
state. While electrical measurements of magnetic resonance have been performed on single spins [iH 
IT^| . the statistical nature of these experiments based on random telegraph noise measurements has 
impeded the readout of single spin states. In this letter, we demonstrate the measurement of the spin 
state of 31 P donor electrons in silicon and the observation of Rabi flops by purely electric means, 
accomplished by coherent manipulation of spin-dependent charge carrier recombination between 
the 31 P donor and paramagnetic localized states at the Si/Si02 interface via pulsed electrically 
detected magnetic resonance. The electron spin information is shown to be coupled through the 
hyperfine interaction with the 31 P nucleus, which demonstrates the feasibility of a recombination- 
based readout of nuclear spins. 



Since the detection of single charges has become tech- 
nically straight forward, it is widely believed 7, 8, 10, 
[HI fl2| that the realization of spin-to-charge transfer is 
the key prerequisite for a successful implementation of 
single spin phosphorus ( 31 P) readout devices, capable of 
determining the actual spin state (spin up | j) or spin 
down I J.)). Different approaches to the electrical spin 
readout of 31 P-donor electron spins have been proposed 
based on spin-dependent transitions between neighbor- 
ing 31 P atoms @ H, EH EI- Since the states involved 
are energetically degenerate, these spin-to-charge trans- 
fer schemes are rather difficult to realize. Alternatively, 
spin-dependent transitions involving dissimilar paramag- 
netic states might be easier to detect as proposed by 
Boehme and Lips 14]. Spin-dependent charge carrier 
transport and recombination are known at least since 
1966 |l5j . when Schmidt and Solomon already observed 
spin-dependent recombination involving 31 P donors in 
silicon. However, it was not demonstrated until 2003 |l6| 
that the much more sensitive electrical detection of spins 
via resonant changes of recombination processes is also 
able to reflect coherent spin motion, which is necessary 
for a readout of the spin quantum state as opposed to a 
mere detection of the presence of spins. 

Figure ^a) illustrates the readout scheme based on 
spin-dependent recombination demonstrated in this pa- 
per. In order to probe 31 P donor electron spins, spin- 
dependent excess charge carrier recombination through 
so-called Pbo centers is used. Pbo centers are trivalent Si 
atoms at the interface between crystalline silicon (c-Si) 
and silicon dioxide (S1O2) that introduce localized, para- 
magnetic states in the Si band gap and dominate elec- 
tron trapping and recombination at the interface [ItI Il8j | . 
If a neutral Pt>o center is located in the vicinity of a 
31 P donor, the electron bound to the 31 P and the electron 



in the Pbo center can form a coupled spin pair |l9l |20( , 
whose wave function |\&) can exist in an arbitrary su- 
perposition of any of its four energy eigenstates. Two 
of these states, \T+) = | TT> and |T_) = |||), where 
the first arrow indicates the state of the 31 P spin and 
the second arrow the state of the Pbo spin, have triplet 
symmetry and are independent of the strength of the 
spin-spin coupling. The remaining two states depend on 
the coupling: For strong spin-spin interaction, they are 
the pure triplet \T ) = ^= (|T4> + lit)) and the singlet 
\S) = -j= (|T-l) — |4T)) states, whereas in the absence of 
spin-spin coupling, they are product states 1 1|) and | J.|) 
which can be represented by equal mixtures of \Tq) and 
\S). For the intermediate coupling regime, the energy 
eigenbase changes continuously between the two limit- 
ing cases. Once a spin pair is generated, it can relax 
by a transfer of the electron from the 31 P donor to the 
Pbo, forming a negatively charged Pbo- This relaxation 
transition depends strongly on the initial symmetry of 
the two spins within the pair. Since the energetically 
lowest doubly occupied charged state of the Pbo has to 
be diamagnetic due to the Pauli exclusion principle, the 
relaxation probability is prop ortional to the singlet con- 
tent I (5'1'Sl/) | 2 of the pair |2fJ. After subsequent capture 
of an excess hole by the P^q, and an excess electron by 
the 31 P + , both centers return to their initial uncharged 
states. This capture of mobile charge carriers e.g. gener- 
ated by illumination renders the spin-dependent process 
detectable by current measurements. 

The incorporation of this readout scheme into an ar- 
ray of 31 P qubits is sketched in Fig. IHb). While cou- 
pling between neighboring qubits is achieved via J gates 
controlling exchange, A gates are used to couple the 
31 P qubits to the Pbo readout centers. Selective read- 
out of a single 31 P in this array is achieved by apply- 
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ing a negative gate potential to all A gates except one 
which is decoupling 31 P and Pbo- At one qubit, a pos- 
itive bias is applied, so that spin-dependent recombina- 
tion exclusively takes place there, monitored by a current 
flowing between the two R contacts, which depends on 
the spin state of the qubit addressed. For an applica- 
tion of spin-dependent recombination to spin readout, 
the 31 P -Pbo mechanism discussed above might be su- 
perior to Kane's original approach 0. Since Pbo states 
have a much higher localization (more than 80% of the 
wave function can be found within a c-Si lattice con- 
stant |21|1 than the 31 P donors [2^, one can expect that 
gate potentials might influence the Pbo wave function 
to a much lesser degree than the donor wave function, 
conceptually simplifying the gate-controlled exchange in- 
teraction. Note that a sufficiently slow application of a 
positive potential on an A gate would not only increase 
the 31 P-Pbo exchange, it would also adiabatically reduce 
the 31 P hyperfine interaction 0,113 ■ Since the nuclear 
spin affects the 31 P-donor electron spin resonance due to 
the hyperfine coupling, the slow application of the posi- 
tive bias would adiabatically encode the 31 P nuclear spin 
state into the permutation symmetry of the 31 P-Pbo elec- 
tron spin pair, while, at the same time, the 31 P nucleus 
would become decoupled. Thus, also a selective and 
relaxation free 31 P nuclear spin-readout determined by 
spin-dependent 31 P-Pbo transitions could be possible via 
selective application of an electric field at an individual 
31 P site in the array in Fig.QJb). 

To demonstrate the feasibility of this readout scheme 
mapping the 31 P donor electron spin state on electrical 
current, pulsed electrically detected magnetic resonance 
experiments (pEDMR) were performed on the test struc- 
ture shown in Fig. ^c). A resonant microwave pulse 
is used to induce coherent manipulation of the electron 
spins by electron spin resonance (ESR). The resulting 
change of the recombination rate is then detected by a 
transient photoconductivity measurement after the co- 
herent excitation is turned off 20] . For the prepara- 
tion of the initial state of the ensemble of spin pairs, the 
sample is allowed to attain a low temperature T = 5 K 
steady state in the presence of a strong offset photocur- 
rent I = 50 /iA, a constant magnetic field Bq ss 350 mT 
and illumination with a tungsten light source. In this 
way, a high density of \T + ) and |T_) spin pairs is gener- 
ated, as the steady state recombination leads to a deple- 
tion of the short lived pair states \S) and, in the case of 
intermediate or weak coupling, \Tq). pEDMR is induced 
by an intense, coherent ESR microwave pulse generated 
by a Bruker E580 pulse bridge and amplified by a 1 kW 
travelling wave tube amplifier, which imposes a unitary 
transformation on the spin pair eigenstates. The final 
state of this transformation is a coherent non-eigenstatc, 
determined by adjustable pulse parameters such as inten- 
sity, pulse length and microwave frequency. If this state 
has an increased singlet content, the recombination rate 




FIG. 1: (a) Sketch of the readout of P spin quantum states 
via spin-dependent charge carrier recombination at Pbo cen- 
ters at the crystalline Si/Si02 interface. The electrons at the 
31 P donor and the Pbo center form a weakly-coupled spin pair 
|^). Since the negatively charged P^-state is diamagnetic the 
recombination probability is proportional to the singlet con- 
tent |(S|\I/)| 2 of the initial pair. Electron spin resonance of 
either constituent of the pair can therefore influence the re- 
combination of excess charge carriers, (b) Incorporation of the 
31 P-Pbo readout scheme into an array of 31 P qubits coupled by 
J gates. Readout of a specific 31 P qubit is achieved by selec- 
tive application of a positive gate bias to the corresponding A 
gate, (c) Sample structure used in the experiments to demon- 
strate electrical detection of Rabi oscillations of 31 P qubits. 
To preferentially locate 31 P in the vicinity of c-Si/SiCh-states, 
the active sample is a 15 nm thick layer of c-Si, doped with 
10 17 cm -3 31 P atoms, on a 5 /im intrinsic buffer epitaxially 
grown by low-pressure chemical vapor deposition on a 30 Qcm 
Boron-doped Si(OOl) wafer. A high density of Pbo centers is 
achieved by a native oxide. The contacts are a grid structure 
consisting of 10 nm thick chromium and 100 nm thick gold 
films and comprising 112 digits of 2 mm length, 10 /im width 
and 10 fim distance. 



after the pulse is also increased and will only gradually 
return to its steady state, as recombination processes of 
the charge carrier pairs typically take place on a much 
slower timescale than their manipulation by the strong 
microwave pulse. Figure Gla) displays a contour plot 
of the change AI of the photocurrent from its steady 
state value as a function of time t after the microwave 
pulse was turned off. The data shown were obtained 
for different magnetic fields Bq and represent the mea- 
sured data after subtraction of non-resonant microwave 
artifacts. The contour plot displays three resonances at 
B = 346.2 mT, B Q = 347.1 mT, and B Q = 350.3 mT, 
also visible in Fig.^b) which shows a time slice through 
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FIG. 2: Contour plot (a) of the transient current change after 
a 480 ns long microwave pulse with a power Pmw = 1 W as a 
function of the magnetic field Bo and after subtraction of mi- 
crowave artifacts. To protect the amplifier against overload, 
the detection system is activated only after 3 [is. Every time 
slice has the shape of an ESR spectrum as exemplarily shown 
for t = 15.5 /is (b). The spectrum is a superposition of the 
hyperfine-split resonance of the 31 P donors and the resonance 
of the Pbo centers. 

the contour plot for t = 15.5 fis. The two features in- 
dicated by dashed lines have equal amplitudes, are sep- 
arated by 4.2 mT and correspond to a central g-value 
of 1.9985. This is the characteristic ESR hyperfine sig- 
nature of the 31 P donor electron in silicon OHil- The 
broader peak at B = 347.1 mT is an unresolved superpo- 
sition of signal contributions from Pbo centers that have 
different orientations with respect to the Si (100) surface, 
which results in resonances at g — 2.0039 and g = 2.0081 
for the orientation of the sample with respect to the mag- 
netic field orientation indicated in Fig. \£\c) 0, ■ The 
current transients at the three resonances exhibit iden- 
tical dynamics: An initial strong photocurrent decrease 
caused by the higher singlet content followed by a slower 
increase of the current originating from the quenching 
of the number of triplet states which also have a finite 
recombination probability |2(ij |. We conclude from the 
identical time dependence that the electronic transitions 
which cause the observed signals must belong to the same 
recombination process, involving 31 P donors and Pbo cen- 
ters, as illustrated in Fig. 

The transient response of the photocurrent after the 
microwave pulse is proportional to the relative singlet 
density change induced by the microwave pulse |20j . 
Therefore, a current change, integrated over a certain 
time in a box-car type of experiment for the purpose of 
noise reduction yielding a charge Q, reflects the spin state 
densities directly after the pulse. Since a change of the 
pulse length r causes a different nutation angle for the 
31 P spin on resonance and thus a different singlet content 
for the 31 P-Pbo spin pair after the pulse, a measurement 
of Q(t) must reveal the coherent Rabi nutation of the 
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FIG. 3: Electrical observation of Rabi flops of 31 P elec- 
tron spins, (a) Photoconductivity change integrated from 
ii = 7 fJ,s to t2 — 23 /is after a coherent ESR excitation 
as a function of the applied pulse length r for four different 
microwave powers Pmw- (b) Plot of the Fast-Fourier Trans- 
form of Q(t) for the four data sets shown in (a), normalized 
to the intensity of the frequency peak, (c) Plot of the peak 
frequencies Q from (b) as a function of the microwave field 
amplitude Bi exhibiting the linear relationship expected for 
Rabi flops. 



spins in resonance. 

To observe Rabi flops of the 31 P spins, the transient 
current measurement was repeated on resonance with the 
31 P peak at the magnetic field of Bq = 350.3 mT un- 
der application of microwave pulses with different length 
r. Figure |^a) displays the measurement of Q(t) ob- 
tained for four different pulse powers Pmw under oth- 
erwise identical experimental conditions as for the ex- 
periments shown in Fig. [21 We clearly see well devel- 
oped oscillations which increase in frequency for increas- 
ing microwave power. This can be seen more quantita- 
tively in Fig.^b), which displays the Fast-Fourier Trans- 
form (FFT) of Q(t). These curves were normalized to 
the maximum of their respective frequency peaks. With 
increasing microwave field strength B\ oc \JPmw , we 
observe a clear shift of the characteristic oscillation fre- 
quency to higher values. Fig.^Jc) shows that the peak 
frequency increases linearly with the excitation field B\ , 
as expected for Rabi oscillations. For the integration 
times used in Fig. [31 a maximum of the recombining 
charge Q(t) corresponds to a maximum of singlet states 
after the pulse. Since at low temperatures the spins of 
31 P and Pbo are in the | [[) triplet state before the pulse, 
the maximum of Q indicates | j) of 31 P after the pulse, 
a minimum corresponds to ||). Rabi echo experiments 
(not shown) reveal that the fast decay of the Rabi os- 
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FIG. 4: Contour plot of the Fast-Fourier Transform of Q{j) as 
a function of the applied magnetic field Bo- The data confirm 
the predictions for the off-resonance Rabi frequency of a 
weakly-coupled spin-1/2 system with g = 1.9985 indicated by 
the white solid line. 

cillations in Fig. |3fa) is not due to incoherence but is 
caused by an inhomogeneous B\ field which leads to a 
fast coherent dephasing of the spin ensemble. 

On resonance, £1 = jBx should hold for isolated spin 
1/2 states, where 7 = gun/h is the gyromagnetic ratio 
and hb the Bohr magneton. From the linear fit of the 
data displayed in Fig. Etc), we obtain a proportionality 
factor of 7 oxp = (1.03 ± 0.08)7. Since the ratio fl/jBi 
is greater than 1 for strongly coupled systems |2(|, the 
experimentally observed fl/^fBi w 1 indicates a weakly 
coupled spin pair. 

Off resonance, the relation between the observed os- 
cillation frequency of Q(t) and the applied B\ field 
is described by the Rabi nutation frequency formula 
f2 = $/J yB] ) 2 + (ui — lot) 2 in the limit of uncoupled 
spins |20ll27j]. Here, lo is the angular frequency of the mi- 
crowave radiation and ojl = jBq the Larmor frequency 
of the spin. To test the Rabi frequency formula for spin- 
dependent recombination off resonance, Q was measured 
with a fixed B\ and varied Bq (or Larmor frequency lol) 
near the high-field 31 P resonance. Figure 0] shows the 
results of these measurements in a contour plot of the 
FFT of Q(t) as a funtion of Bo- A symmetric increase 
of fl with increasing distance from the resonance peak is 
observed in perfect agreement with the predictions for 
at g = 1.9985, indicated by the white line. Therefore, the 
oscillations of Q(t) shown in Figs. Eland 01 fully meet the 
predictions of Rabi oscillations of weakly coupled spin 
1/2 states. 

The experiments presented here have been performed 
on ensembles of 31 P donors. The charge noise achieved 
so far is about 10 4 e as shown in Fig. 13 Since the flip of a 
single spin approximately leads to one elementary charge 
detected in the box-car pEDMR experiment [2(| , the de- 
tection limit of pEDMR demonstrated here is around 11 



orders of magnitude lower than for conventional nuclear 
magnetic resonance and about 7 orders of magnitude 
lower than for ESR 28]. The noise most likely is deter- 
mined by the shot noise of shunt currents due to excess 
charge carriers which diffuse into the c-Si bulk. Suppres- 
sion of these currents via improved sample geometry e.g. 
using a buried oxide and lateral structuring for confine- 
ment should reduce the noise level, resulting in an even 
higher sensitivity. In fact, cw-EDMR was recently used 
to successfully detect as few as 100 31 P atoms implanted 
into intrinsic c-Si [2^. 

The experimental results shown are an important step 
towards the realization of an electrical spin readout of 
donors. The experiments demonstrate (i) a spin-to- 
charge transfer mechanism required for 31 P spin readout 
via 31 P-Pbo pairs and show that (ii) this spin pair is suf- 
ficiently weakly coupled so that the 31 P hypcrfinc inter- 
action can be observed, ultimately allowing the electrical 
detection of the nuclear spin state [3(J- The advantage 
of the concept is that it utilizes the native paramagnetic 
state at the c-Si/Si02 interface. The average distance 
of these Pbo centers can be adjusted between 10 nm and 
100 nm by varying oxidation conditions. Given the local- 
ization length of the 31 P donor [22J, it is therefore possi- 
ble to prepare an interface where each 31 P qubit interacts 
with one Pbo center. 

This work was funded by Deutsche Forschungsgemein- 
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grown by G. Vogg and F. Bensch at Fraunhofer IZM in 
Munich. 
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